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Abstract 

We investigate thermodynamic phases, including the phase of coexistence of 
superconductivity and ferromagnetism, the possible phase transitions of first and 
second order, and the shape of the phase diagram in mean-field approximation 
for a phenomenological model of spin-triplet ferromagnetic superconductors. The 
results are discussed in view of application to metallic ferromagnets as UGe2, ZrZn2, 
URhGe, and Fe. 

Recently, the coexistence of ferromagnetism and superconductivity was discovered in the 
metallic compounds UGe2 [UGIE] ZrZn2 @J, URhGe [E] and also in Fe jHj in experiments 
at low temperatures and high pressure. In contrast to other superconducting materials, 
in these metals the phase transition temperature to the ferromagnetic state is higher than 
the phase transition temperature to the superconducting state. Moreover, it seems that 
the superconductivity in the metallic compounds mentioned above always coexists with 
the ferromagnetic order and is enhanced by the latter. In these systems the superconduc- 
tivity seems to arise from the same electrons that create the band magnetism, and is most 
naturally understood as a triplet rather than spin-singlet pairing phenomenon EJ IH] • 
The same unconventional superconductivity has been suggested [l] as a possible outcome 
of the interpretation of experiments in Fe 0. Note, that both vortex and Meissner su- 
perconductivity phases jE] are found in the high-pressure crystal modification of Fe which 
has a hexagonal close-packed lattice. In this hexagonal lattice the strong ferromagnetism 
of the usual bcc iron crystal probably disappears [TU] . 

Recently, the phenomenological theory that explains coexistence of ferromagnetism and 
unconventional spin-triplet superconductivity of Landau-Ginzburg type was developed 
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[TT1^|. The possible low-order coupling between the superconducting and ferromagnetic 
order parameters is derived on the basis of general symmetry group arguments and several 
important features of the superconducting vortex state in the ferromagnetic phase of 
unconventional ferromagnetic superconductors are established I12j. 

In this letter we shall use the approach presented in Refs. [TT) IT2*] to investigate the con- 
ditions for the occurrence of the Meissner phase and to demonstrate that the presence 
of ferromagnetic order enhances the p-wave superconductivity. For this aim we shall es- 
tablish the phase diagram corresponding to model ferromagnetic superconductors in a 
zero external magnetic field. We shall also show that the phase transition to the super- 
conducting state in ferromagnetic superconductors can be either of first or second order 
depending on the particular substance. We confirm the predictions made in Refs. [THIT2] 
about the symmetry of the ordered phases. 

Our investigation is based on the mean-field approximation as well as on familiar 
results about the possible phases in nonmagnetic superconductors with triplet (p-wave) 
Cooper pairs pHHUHH]. We shall neglect all anisotropics, usually given by the respective 
additional Landau invariants and gradient terms jSHH] in the Ginzburg-Landau free energy 
of unconventional superconductors. The reasons is that the inclusion of crystal anisotropy 
is related with lengthy formulae and a multivariant analysis which will obscure our main 
aims and results. Let us emphasize that the present results should be valid in the same 
or modified form when the crystal anisotropy is properly taken into account. We have to 
mention also that there is a formal similarity between the phase diagram obtained in our 
investigation and the phase diagram of certain improper ferroelectrics [TT)j I17j. 

Following Refs. fUEJ we consider the Ginzburg- Landau free energy F = J d 3 xf(if>, A4), 
where 

/ = ^-{D*^){D^) + a s H 2 + b -\^ + a f M 2 + ^M 4 + i lo M.& x V*) • (1) 

In Eq. (|TJ), Dj = (V — 2ieAj/hc), where Aj (j = 1,2,3) are the components of the 
vector potential A related with the magnetic induction B = V x A , the complex vector 
ip = {ipj} = i>2, ^3) is the superconducting order parameter, corresponding to the 
spin-triplet Cooper pairing and M. = {Aij} is the magnetization. The coupling constant 
70 = 47rJ > is given by the ferromagnetic exchange parameter (J > 0). Coefficients 
a s = ot s {T — T s ) and a/ = a/(T — Tj) are expressed by the positive parameters a s and 
«/ as well as by the superconducting (T s ) and ferromagnetic (T/) critical temperatures 
in the decoupled case, when TW^-interaction is ignored; b > and (3 > as usual. Note, 
that all seven material parameters (a s , a/, T s , T s , b, /3, J) depend on material properties, 
the temperature T and additional intensive thermodynamic parameters as pressure P. 

We assume that the magnetization M. is uniform, which is a reliable assumption outside 
a quite close vicinity of the magnetic phase transition whereas we keep the spatial (x—) 
dependence of if). The reason is that the relevant dependence of if) on x is generated by 
the diamagnetic effects arising from the presence of Ai and the external magnetic field 
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H 112] rather than from fluctuations of ip (this effect is extremely small and can be 
safely ignored). Note, that the first term in (0) still persists for H = because of the 
diamagnetic effect created by of magnetization M. = B/Att > 0. As we shall investigate 
the conditions for the occurrence of the Meissner phase where ip is uniform, the spatial 
dependence of ip an d, hence, the first term in r.h.s. of ((TJ will be neglected (see also a 
brief discussion at the end of this paper). 

In case of a strong easy axis type of magnetic anisotropy, as is in UGe2 jl, the overall 
complexity of mean-field analysis of the free energy can be avoided by performing 
an Ising-like description: M. = (0,0, M), where M. = ±|.M| is the magnetization along 
the "z-axis." Further, because of the equivalence of the two physical states (±.M) the 
thermodynamic analysis can be performed within the "gauge" M. > 0. But this stage 
of consideration can also be achieved without the help of crystal anisotropy arguments. 
When the magnetic order has a continuous symmetry one may take advantage of the 
symmetry of model and avoid the consideration of equivalent thermodynamic states 
that occur as a result of the respective symmetry breaking at the phase transition point 
but have no effect on thermodynamics of the system. In the isotropic system one may 
again choose a gauge, in which the magnetization vector has the same direction as z-axis 
(\M\ = M-z = M.) and this will not influence the generality of thermodynamic analysis. 

With the help of convenient notations, ipj = ft^ipj, V 9 ? = 07 ex P(#j)> M = /3 1//4 .A/1, 7 = 
7o/(6 2 /?) 1 / 4 , r = a s /Vb, t = af/yf]3, and neglecting the first term in ((TJ) the free energy 
becomes 

/ = r\<j)\ 2 + -|0| 4 + tM 2 + -M 4 + 2 7 M0 1 2 sin# , (2) 
where 2 = (0 2 + 2 + 2 ), and 6 = {6 2 - 9 X ). 

The possible (stable, metastable and unstable) phases are given in Table 1 together with 
the respective existence and stability conditions. The stability conditions define the do- 
main of the phase diagram where the respective phase is either stable or metastable jj.13, . 
The normal (disordered) phase, denoted in Table 1 by iV always exists (for all tempera- 
tures T > 0) but is stable for t > 0, r > 0. The superconductivity phase denoted in Table 
1 by SCI is unstable. The same is valid for the phase of coexistence of ferromagnetism 
and superconductivity denoted in Table 1 by C02. The N-phase, the ferromagnetic phase 
(FM), the superconducting phases (SC1-3) and two of the phases of coexistence (COl-3) 
are generic phases because they appear also in the decoupled case (7 = 0). When the 
M^-coupling is not present, the phases SC1-3 are identical and represented by the order 
parameter ip where the components ipj participate on equal footing. The asterisk attached 
to the stability condition of "the second superconductivity phase" ,(SC2), indicates that 
our analysis is insufficient to determine whether this phase corresponds to a minimum 
of the free energy. As we shall see later the phase SC2, as well as the other two purely 
superconducting phases and the coexistence phase COl, have no chance to become stable 
for 7 7^ 0. This is so, because the non-generic phase of coexistence of superconductivity 
and ferromagnetism (FS in Table 1), which does not exist for 7 = is stable and has a 
lower free energy in their domain of stability. 
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TABLE 1. Phases and their existence and stability properties (k = 0, ±1, ...). 
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We have outlined the domain in the (t, r) plane where the FS phase exists and is a 
minimum of the free energy. For r < the third-degree algebraic equation for M (see Table 
1) and the existence and stability conditions are satisfied for any M > provided t > 7 2 . 
For t < 7 2 the condition M > M have to be fulfilled, here the value M = (— 1 + 7 2 /2) 1 / 2 
of M is obtained from r(Mo) = 0. Thus for r = the N-phase is stable for t > 7 2 /2, on 
the other hand FS is stable for t < 7 2 /2. For r > 0, the requirement for the stability of 
FS leads to the inequalities 



(3) 



max [ -,M m ) < M < M , 



where M m = (M /y/3) and M should be the positive solution of the third-degree equation 
of state from Table 1; M m > gives a maximum of the function r(M). 

The further analysis leads to the existence and stability domain of FS below the line AB 
given by circles (see Fig. 1). In Fig. 1 the curve of circles starts from the point A with 
coordinates (7 2 /2, 0) and touches two other (solid and dotted) curves at the point B with 
coordinates (— 7 2 /4, 7 2 /2). Line of circles represents the function r(M m ) = r rn (t) where 



r m (t) 



4 (l 2 
3^7 V 2 



-t 
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(4) 



Dotted line is given by r e (t) = 7y The inequality r < r m (t) is a condition for the 
stability of FS, whereas the inequality r < r e (t) for (—t) > 7 2 /4 is a condition for the 
existence of FS as a solution of the respective equation of state. This existence condition 
for FS has been obtained from 7M > r (see Table 1). 

In the region on the left of the point B in Fig. 1, the FS phase satisfies the existence 
condition 7M > r only below the dotted line. In the domain confined between the lines 
of circles and the dotted line on the left of the point B the stability condition for FS is 
satisfied but the existence condition is broken. The inequality r > r e (t) is the stability 
condition of FM for < (—t) < ^i 2 /A. For (—t) > 7 2 /4 the FM phase is stable for all 
r > r e (t). 
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Figure 1: The phase diagram in the plane (t, r) with two tricritical points (A and B) and 
a triple point C; 7 = 1.2. 



In the region confined by the line of circles AB, the dotted line for < (—t) < 7 2 /4, and 
the t— axis, the phases N, FS and FM have an overlap of stability domains. The same 
problem exists for FS and the SC phases in the second quadrant and for the phases FS 
and COl in the third quadrant of the plane (t, r). The comparison of the respective 
free energies for r < shows that the stable phase is FS whereas the other phases are 
metastable within their domains of stability. 

The part of the t-axis given by r = and t > 7 2 /2 is a phase transition line of second 
order which describes the N-FS transition. The same transition for < t < 7 2 /2 is 
represented by the solid line AC which is the equilibrium transition line of a first order 
phase transition. This equilibrium transition curve is given by the function 

r eq (t) = ][3 1 -( 1 2 + lQt)) 1 / 2 ]M eq (t), (5) 



where 



1 



1/2 



M eg (t) = [ 7 2 - 8t + 7 ( 7 2 + 16t) 1/2 ] ~ (6) 

is the equilibrium value (jump) of the magnetization. The order of the N-FS transition 
changes at the tricritical point A. 

The domain above the solid line AC and below the line of circles for t > is the region 
of a possible overheating of FS. The domain of overcooling of the N-phase is confined by 
the solid line AC and the axes (t > 0, r > 0). At the triple point C with coordinates [0, 



r eg(0) = 7 2 /4] the phases N, FM, and FS coexist. For t < the straight line 

K q (t) = ^ + \t\, - 7 2 /4<t<0, (7) 

describes the extension of the equilibrium phase transition line of the N-FS first order 
transition to negative values of t. For t < (— 7 2 /4) the equilibrium phase transition FM- 
FS is of second order and is given by the dotted line on the left of the point B (the 
second tricritical point in this phase diagram). Along the first order transition line r* (t) 
given by (jHJ) the equilibrium value of M is M eq = 7/2, which implies an equilibrium order 
parameter jump at the FM-FS transition equal to (7/2 — yjtj). On the dotted line of 
the second order FM-FS transition the equilibrium value of M is equal to that of the FM 
phase (M eq = \/Jt|). Note, that the FS phase does not exists below T s and this seems to 
be a disadvantage of the model (0). 

The equilibrium phase transition lines of the FM-FS and N-FS transition lines in Fig. 1 
can be expressed by the respective equilibrium phase transition temperatures T eq defined 
by the equations r e = r(T eq ), r eq = r(T eq ), r* eq = r(T eq ), and with the help of the relation 
M eq = M{T eq ). This leads to some limitations on the possible variations of the parameters 
of the theory. For example, the critical temperature (T eq = T c ) of the FM-FS transition 
of second order (7 2 /4 < —t) is obtained in the form T c = (T s + AirJJ\A/a s ), or, using 

T C = T S + 

2 

where Tf > T s , and T* = (47[J) 2 a.f/a 2 s f3 is a characteristic temperature of the model (JTJ). 
The investigation of the conditions for the validity of (jHJ) leads to the conclusion that the 
FM-FS continuous phase transition (at 7 2 < —t) is possible only if the following condition 
is satisfied: 

T f - T s > = (q + ^)T* , (9) 

where ? = j3a^/4ba'j. This means that the second order FM-FS transition should disap- 
pear for sufficiently large M-^-coupling. Such a condition does not exist for the first order 
transitions FM-FS and N-FS. 

Taking into account the first term in the free energy (0) should lead to a depression of 
the equilibrium transition temperature. As the magnetization increases with the decrease 
of the temperature, the vortex state should occur at temperatures which are lower than 
the equilibrium temperature T eq of the homogeneous (Meissner) state. For example, the 
critical temperature (T c ) corresponding to the inhomogeneous (vortex) phase of FS-type 
has been evaluated [TJ] to be lower than the critical temperature (JHJ): (T c — T c ) = 
47thb-M/c( s (fis — \e\h/2mc - Bohr magneton). For J ^> /ig, we have T c ps T c . Finally, 
let us emphasize that a more reliable description of these phenomena, in particular, of 
the thermodynamic behaviour of the FS phase at relatively large values of M. can be 
performed if an additional term of type .M 2 |?/>| 2 is included in the model (JT|). 



(rp* \ 2 
T J +T*(T f -T s 
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